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Thomas W. Orange 
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SUMMARY 

A method Is presented for estimating the crack-extension resistance curve 
(R-curve) from residual-strength (maximum load against original crack length) 
data for precracked fracture specimens. The method allows additional Informa- 
tion to be Inferred from simple test results, and that Information can be used 
to estimate the failure loads of more complicated structures of the same 
material and thickness. 

The paper first reviews the fundamentals of the R-curve concept. Then the 
analytical basis for the estimation method Is presented. The estimation method 
has been verified In two ways. Data from the literature (Involving several 
materials and different types of specimens) were used to show that the esti- 
mated R-curve Is In good agreement with the measured R-curve. A recent pre- 
dictive blind round-robin program offered a more crucial test. When the 
actual failure loads were disclosed, the predictions were found to be In good 
agreement. 


INTRODUCTION 

The crack-extension resistance curve, or R-curve, Is one of the most pow- 
"'erful concepts available to the fracture analyst. In this paper, a useful 
extension to that concept will be described. The R-curve can be used to pre- 
dict failure loads for any Initial crack size In any specimen or structural 
configuration (of the same material and thickness and In the same environment) 
for which a stress Intensity analysis Is available. Now It Is possible to 
reverse the process. That Is, If the residual strength curve (fracture stress 
against original crack size) Is known, the R-curve can be determined. 

First the fundamentals of the R-curve concept will be reviewed. Although 
It has been presented earlier, the derivation of the R-curve estimation method 
will also be reviewed. Then three applications of the estimation method will 
be described. The method was used to show that several semlempirical fracture 
analyses are each equivalent to a particular R-curve formulation. The R-curve 
was estimated to a useful degree from published residual strength data for a 
variety of specimen configurations. Finally, the estimation method was tested 
In an ASTM predictive blind round-robin program. 


R-CURVE CONCEPT 

Four papers are significant In the history of the R-curve concept. Irwin 
(ref. 1) first Introduced the concept In 1954 and corrected It In 1959 
(ref. 2). In 1961 , Krafft (ref. 3) extended the concept by postulating that 
(a) for a given material and thickness, crack propagation resistance depends 
only on crack extension, and (b) an effective crack extension can be Inferred 
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from compliance measurements. He stated that "...the information necessary to 
predict fracture Instability oyer a wide range of dimensions Is not a single 
Gc value but the shape of the entire resistance curve." In 1968, Clausing 
(ref. 4) showed that "...the complete 6 r curve can be determined In one 
stable specimen by measuring load and crack length as the crack propagates," 

The R-curve concept Is Illustrated schematically In figure 1(a) for an 
Infinite body containing a crack whose original length Is ZBq. The strain 
energy release rate, G^, Is given by 



where a Is the applied stress and E* Is the effective modulus. G^ 
represents the driving force (per unit thickness) tending to cause crack prop- 
agation. The material's resistance to crack propagation, Gr, Is a function 
of crack extension A. As stress normal to the crack Is applied and Increased 
to 90 percent of the critical stress In figure 1(a), the crack must extend 
only a small distance to develop a large resistance. At this point the crack- 
extension resistance equals the driving force and the crack Is stable. As the 
stress Is Increased, progressively larger amounts of crack extension are 
required to resist the crack driving force. Finally at the critical stress oq 
the driving-force curve and the R-curve are tangent. Beyond the point of tan- 
gency the driving force Increases faster with crack length than does the mate- 
rial's resistance. This Instability condition represents the failure of the 
body. The point of tangency defines the fracture toughness Gj. and the 
critical crack length 2aQ. Since the driving force curve for an Infinite 
body Is a straight line, both the fracture toughness and the amount of crack 
extension at Instability Increase with Increasing original crack length. If 
the R-curve exhibits a plateau, Gj. and Aj. may asymptotically approach 
limit values. 

In simple finite bodies and test specimens the presence of stress-free 
boundaries results In an additional Increase In the crack driving force as the 
crack extends toward a boundary. Thus the slope of the driving-force curve 
Increases continuously with Increasing crack length. The Instability condition 
for a typical finite-width specimen Is shown In figure 1(b). As the Initial 
crack length Is Increased from zero, both Gj. and Aj- Increase at first, 
but reach maximum values that depend on the specimen width and on the forms of 
both the driving force curve and the R-curve. As the Initial crack length Is 
Increased still further, both Gc ^nd Aj. begin to decrease. 

Instability predictions can be done graphically or analytically. Graph- 
ical methods will not be treated here, but the analytical method Is outlined 
In the following section. 


METHOD OF ANALYSIS 

Conventional Instability calculations require differentiable mathematical 
expressions for both the crack driving force curve and the R-curve. The former 
Is obtained from elastic fracture mechanics theory, the latter by curve-fitting 
to experimental data. For Instability to occur, the magnitudes and the slopes 
of both curves; must be equal . This requi res the simultaneous solution of a 
pair of equations. As shown In references 5 and 6, following some algebraic 
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manipulation, the Instability condition can be written In a particularly useful 
form as 


0 = 


9(‘c> ‘o * ‘e 

g'(»,) - 1 * 2., 


where g(A) = E'Gr 

g'(A) = E' dGp/dA 
a = (a/Y)(dY/da) 


( 1 ) 


and A Is the effective crack extension, a^ Is the Initial crack length, 

Y Is the stress Intensity calibration factor, and the subscript (c) means 
"evaluated at the Instability point." This form Is useful for computing since 
the first term Includes only material properties and the second term Includes 
only geometrical parameters. Now If the functions g(A) and g'(A) and the 
appropriate equation for a are substituted Into this equation then, for 
prescribed values of ag and W, Ag is the least positive root of that equa- 
tion. This root can be found by any of several numerical methods. Finally, 
the fracture stress Is calculated from the crack driving force expres- 
sion. By repeating the calculation for relative crack lengths from near-zero 
to near-unity, the complete residual strength curve may be developed 
poInt-by-poInt. 


Conversely, If an equation for the residual strength curve Is available. 

It Is possible to derive an expression for the_ corresponding R-curve. To do 
this we must differentiate one of the Instability equations with respect to the 
critical crack extension Ag. That operation Is described In detail In 
references 5 and 6, so only the results will be given here. The mathematics 
are greatly simplified If we reduce the problem to one of a single Independent- 
variable. This Is done by prescribing the manner In which ag and W may 
vary. Three cases will be considered. 


for W = constant 


Case I : W = Constant 

Assume that there Is a function f such that we can define 

f(ag) = 

f'(ag) = df(ag)/dag 

Following differentiation and algebraic manipulation we have 

0 = (1 + 2ag)f(ag) + (ag + A(.)f{ag) 

The use of this equation will be explained shortly. 


(2a) 
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, : Case ll:' ao/W = Constant 

Now assume that there 1s a function h such that 

h(ao) = o? 
h'(ao) = dh(ao)/dao 
and corresponding to equation (2a) we have 

A 


for ap/W = constant 


= (’ - 


)> (a^ + A^) h'(a^j) 


(2b) 


Case III: ao = Constant 


This case Is 
Assume that there 


of limited usefulness but Is Included for completeness. 
Is a function j such that 


J(W) = 

J'(W) E dJ(W)/dW 


for 


ao = constant 


and corresponding to equations (2a) and (2b) we have 




ilM 

J(W) 


(2c) 


The R-curve Is developed as follows. First, at a given point on the residual 
strength curve, the appropriate function (f, h, or j) and Its derivative are 
computed and substituted In the appropriate one of equations (2a, b, and c), 
which must then be solved for A(.. A numerical solution Is usually 
required. Then A^ and <j^ are used to calculate Gq. The 
resulting pair (G(., A(.) define a point on the R-curve. The process Is 
repeated at additional points on the R-curve. If desired, the R-curve can be 
expressed In terms of KR = (E'Gr)^/ 2, Also, a suitable function can be 
fit to describe Gr or Kr as an explicit function of A. 


APPLICATIONS 
Equivalent R-Curves 

In reference 5 this method was used to examine the relations between sev- 
eral semlempirical fracture analyses (SEFA) and the R-curve concept. That 
study explained why a SEFA might yield good results with one set of data and 
poor results with another. It was shown that for each SEFA there Is an equiv- 
alent R-curve whose magnitude and shape are determined by the SEFA formulation 
and Its empirical parameters. That R-curve Is equivalent In that It predicts 
exactly the same relation between fracture stress and Initial crack size 
(residual strength) as the SEFA. 


This was done by rewriting the SEFA formulation (If necessary), then sub- 
stituting Into equation (2a), which can then be solved to give an explicit 
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formulation for the effective R-curve. The equivalent R-curves were found to 
differ markedly from one SEFA to another. Figure 2 (from ref. 5) shows (In 
dimensionless form) equivalent R-curves for four SEFA, (refs. 7 to 10). The 
Newman equivalent R-curve Is asymptotic to unity In figure 5(c). The Bockrath 
curve Is exponential and thus has no asymptote. The shapes of the equivalent . 
R-curves are Intrinsic to the specific SEFA formulation. 

It was concluded that If, for a given set of data, a SEFA correlates 
residual strength closely. Its equivalent R-curve will closely approximate the 
effective R-curve :of the material. It was further concluded that, of the five 
SEFA examined, Newman's (ref. 9) appears to be the most generally useful. In 
Ref. 2, Newman's SEFA was used to generate equivalent R-curves, which were 
then used to predict the load at 5 percent secant offset for center-crack 
(fig. 3(a)) and compact (fig. 3(b)) specimens (figures from ref. 5, data from 
the literature) using only the Initial crack length and maximum load data. 


Estimation of R-Curve from Residual Strength Data 

In reference 6 this method was extended to allow the estimation of the 
R-curve from residual strength data, without recourse to any SEFA. Here the 
functions f, h, and J In equation (2) are obtained by numerical approxima- 
tion and differentiation of residual strength data. This procedure Is not 
always straightforward, but some comments and hints are given In reference 6. 

Since residual strength tests are almost always done on center-crack spec- 
imens, the examples In reference 6 are limited to that geometry. Figures 4(a) 
and (b) show R-curves estimated from residual strength data for specimens of' 
2219^'T87 aluminum and AM355 steel respectively (figures from ref. 6, data from 
the literature). In figure 4(a), essentially the same R-curve Is obtained from 
specimens of two different widths. In figure 4(b), essentially the same 
R-curve Is obtained from a series of specimens having constant crack length as 
from a series with constant ratio of crack length to specimen width. 

In addition to the limitations Inherent In numerical analysis, this method 
has the same limitation as a standard R-curve test run under load control. 

That Is, for a given specimen the R-curve can only be measured (or estimated) 
to the point at which the crack becomes unstable. Tests run under displacement 
control allow a much larger portion of the R-curve to be measured. 


ASTM Predictive Analytical Round-Robin Program 

Arecent predictive blind round-robin program (ref. 11) offered a more 
crucial test. Complete test data (Including measured R-curves) for compact 
tension specimens of three materials were supplied to all participants. How- 
ever, the R-curves were estimated using only the failure loads and original 
crack lengths. Failure loads were then predicted for additional compact spec- 
imens with different dimensions, for center-crack specimens , and for special 
three-hole specimens which simulate the behavior of a cracked stiffened panel. 
Details of the author's participation are given In Appendix K of reference 11. 
When the actual failure loads were disclosed, the predictions were found to be 
In good agreement. 
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Figures 5(a) to (c) show the measured R-curves supplied to the partici- 
pants and the R-curves estimated using reference 6 for the three test materi- 
als, 7075-T651 and 2024-T351 aluminum alloys and AISI 304 /stainless steel, 
respectively. The bold curves were fitted to the data by the round-robin 
organizers. The fine solid curves represent functions fitted to points on the 
estimated R-curve, the fine dashed curves represent extrapolation of those 
functions beyond the point of Instability. Within the range of applicability, 
the estimated R-curves for the aluminum alloys are In good agreement with the 
measured R-curves. For the very ductile stainless steel, the comparison Is not 
as good. The curves have the same general magnitude but completely different 
shapes. As discussed In the body of reference 11, the 304 specimens exhibited 
very large deformations along the crack line during fracture. Clearly the 
limits of linear elastic fracture mechanics were severely violated, and It Is 
surprising that any R-curve method worked at all. 

Because the purpose of my participation was to test the applicability of 
the estimation method, the limits of the method were strictly adhered to. That 
Is, no attempt was made to predict failure loads by other methods such as a 
limit load analysis, even when It was obvious that limit load would be the 
failure load. For cases where predictions were made, these predictions were 
quite good. As shown In table XVI of reference 11 (as Participant number 11), 
the standard error for all predictions made was less than 0.1. 


SUMMARY 

This method of R-curve estimation appears to be quite useful. It allowed 
an explanation and Interpretation of several semlempirical fracture analyses 
based on the more modern R-curve concept. But more Important. Is the ability 
to estimate the R-curve from residual strength data. This means that addi- 
tional Information can now be Inferred from the simple residual strength tests 
which were run many years ago. Such data may be found In the literature or 
may exist In company files. Now the R-curve can be estimated from that data 
and used to predict the failure loads of more complicated structures (such as 
reinforced panels or cracked holes) of the same material and thickness. 
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Figure 1. - Schematic representation of R-curve instability concept 
for infinite and finite bodies. 



(b) Finite body. 
Figure 1. - Concluded. 


DIMENSIONLESS CRACK EXTENSION RESISTANCE, E'KJ^Gn DIMENSIONLESS CRACK EXTENSION RESISTANCE, 


3.375 



DIMENSIONLESS CRACK EXTENSION, (Oy/Kfl'^A 
(c) Newman's analysis (ref. 9). 


DIMENSIONLESS CRACK EXTENSION. (Ou/Kyc)*-^ A 
(d) Bockrath's analysis (ref. 10). 


Figure 2. - Dimensionless R-curves equivalent to various semiempirical fracture analyses for the case of a crack in an infinite 


liinra 




STRESS AT 5 PERCENT SECANT OFFSET, 



INITIAL CRACK HALF-LENGTH, a^, cm 



0.1 2 3 4 5.6 7 8 

INITIAL CRACK HALF-LENGTH, a^, in. 

(a) Stress at 5 percent secant offset for 7075-T7351 aluminum alloy 
center-crack specimens. 

Figure 3, - 5% secant-offset loads for center-crack and compact specimens. 
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Effective crack extension, A, in. 

(a) R-curve for 2219-T87 aluminum alloy. 

Figure 4. - R-curves estimated from residual strength data. 
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Effective crack extension, A, in. 
(b) R-curve for AM355 steel. 
Figure 4. - Concluded, 
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